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Abstract 

In this paper we present the notions of trail (pseudo-)division, 
generahzed subresultants and generalized subresultant algorithm. 

1 Trail pseudo-division 

We will work in some polynomial ring K[x]. First, we define full polynomial 
as one with non-zero trail coefficient. All through the paper we will deal 
with full polynomials. Let we have two full polynomials / and g. Now we 
describe the process of trail (pseudo-)division. The usual (pseudo-) division 
can be illustrated in the following scheme, where the coefficients of both 
polynomials are written from the left to the right by decreasing powers of 
X and the second, third, etc. lines subtracts one after one from the first 
after some multiplications ("*" mean some coefficient, no other comments 
are needed there): 
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We introduce the trail (psudo-) division, which can be analogously illustrated: 

:^ :^ :^ :^ :^ 
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* * * * x'^g 

Here the eliminations perform in the trail part of "bigger" polynomial / by 
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the trail part of the "smaller" one g. In each step one get polynomial from 
the ideal {f,g) (it has zero in some lower terms). After removing of some 
maximal possible power of x {x ^ {f,g), if deg^{gcd{f , g)) > as / and g 
are full and if gcd(/, g) = 1 then resulting polynomial is obviously belongs 
to {f,g)) one can get the polynomial of degree less than g has. 

This way of division is usefuU in the case of pseudo-division. In the usual 
pseudo-division the first polynomial / is multiplied by some degree of the 
leading coefficient of the second one g. In the trail division it is multiplied 
by some degree of the trail coefficient of g. If it is "less" in some sense 
than the leading coefficient then the resulting trail pseudo-remainder will 
have "smaller" coefficients than the usual one. Analogously in the case of 
division (not pseudo) the dividing of g can be performed by " smaller" term. 
The tprem{f, g) will denote the trail pseudo remainder of / and g. The 
remark: in general we can change the "place" of elimination (e.g first vanish 
the leading coefficient, then the trail, again the trail, etc.). 

If h is full polynomial than h* will denote the reverted polynomial (e.g. 
{5x'^ + Ax^ + 3x'^ + 2x + l)* = x^ + 2x^ + + 4x + 5). The following formula 
is valid up to multiplying by some power of x : tprem(f, g) = (prem(f*, g*))*. 
To prove this formula one need to place the mirror near the scheme for trail 
pseudo-remainder. Then in the mirror one will see the process of finding 
usual pseudo-remainder of reverted polynomials. 

Our next goal is to develop the algorithm analogous to the subresul- 
tant algorithm for gcd computation with using trail pseudo-remainders. 
For this purpose we fix here the generalized algorithm for pseudo-remainders 
genPRem: this algorithm gets as input two full polynomials / and g, deg^g < 
deg^f. In output it produces the polynomial r together with the following 
six values: r, 6, A, g, g, w. r is the full part of trail or usual pseudo re- 
mainder depending on which way is better (usaul pseudo-remainder algo- 
rithm doesn't exclude "superflous" powers of x); 6 = deg^f — deg^g; A = 
trailDeg2,(prem(f, g)) if usual pseudo-division is used and trailDeg2.(prem(f*, g*)) 
if trail pseudo- division was performed; g and g are \cxg and tCxg or tCxg and 
Ic^g depending on the way of division: first pair in the usual case and the 
second in the trail one; w is marker of kind of division: lead or trail. For- 
mally we can write it in the following way: 

Algorithm genPRem 

Input: u,v are full polynomials, deg^(u) > deg^(f) 
Output: the generalized pseudo- remainder 
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if relativeSize(lCx(u)) < relativeSize(tCx(u)) then 
w :— prem(u, v); 

return('u; / x^"^^^^^^^^ ,deg^ (u) — deg^. (v) jtrailDeg^. (w) ,lcx (u) ,tCx (u) ,lead) 
else 

w :— prem(u*, v*); 

return( (w/x*''^™^^^ ('"^ ) * ,degj. (li) - degj. (t;) jtrailDegj. (w) ,tCj; (li) ,lca; (li) ^ 

fi; 

Here relativeSize is a integer characteristic of some term which says how 
big it is. For example, the amount of memory which takes the term can be 
used. 



Generalized subresultant algortihm and gen- 
eralized subresultants 



Let f,ghe the initial full polynomials, deg^f > deg^g. Let ui = ui — 
Ui = f, U2 = U2 = U2 = g, 113,114, ... be the sequence of generalized re- 
mainder: Ui = genPRem(ui_2, iii_i). Of course, the elements of this se- 
quence contains removable factors, we need this sequence just to define 
the sequence Sf. we denote 5i — deg^Ui+i — degj.Ui, — YA=m^u — 
X — value of genPRem(ui_i, ili). As in the subresultant algorithm we will 
investigate the determinants of matrices which consist of coefficients of poly- 
nomials a;"/, x'^g: let 
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fixed a cols 



fixed 2(fe+l)+i5i-a-l cols 



where k < deg^g. We denote by {ui, M2)"' the polynomial whose coefficients 
are obtained by fixing some a < rows(Mj) columns in the left part of Mj, 
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rows(Mj) — a — 1 columns in the right side of and calculating the deter- 
minant which consist of fixed columns and one non-fixed column. We are 
not specifying what is a in our consideratuions, we know that there is some 
a. ti2)'^2 will be denoted as Uk+2- Our goal is to express such polynomial 
via taking generalized pseudo-remainders. Most of equations bellow will be 
true up to the sign - the sign is not important in our considerations and 
it's checking is redundant. The following relation will be usefuU for us: it 
describes what is happened when we perform the generalized pseudo-division 
in the matrix: 

From this formula we see, for example, that U4 = (ui,U2)^^ = %7 — j7^M4, 

92 92(52 ) 2 

where M4 = genPRem(ii2, U3), M3 = genPRem(ui, U2), i.e. we know what the 
expression can be removed from M4. Let Ui denote genPRem(iii_2, Ui-i). We 
want to determine how Ui linked with -Uj. Let us fix the number k. Then we 
can write down the following sequence of equations: 

Now we proceed the same transformations with k instead of — 1 and 
simultaneously we will express Gf"*"^ via Gf"*"^ using the ((H): 



U.+2 = (ni,n2)^^^ = |^^Gr^(^2,t.3)^^ 

.U9;^ 



'^k+l[^k,Uk+l) —Ljk+2^k+2- 



^9^-1 



Hence 
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9I' 9u\Y{G]g]L-r"GlX\gt' ) 
Let us denote the expression with product as hk+2'- 

hu^2 = \{G]c^-r'G^^-'. (2) 



hk+2 is "integer" as it is equal to the determinant with "integer" entries: 

= I[G^9f:r'G'^iuk-s,u,.2)'^-^-' = l[G^gfzr'Gt'iuk-2,u,^i)''-^-' 

= I[G]9'-r'G'^9t'-'uk, (3) 

and taking the leading or trail coefficient we get hk+2- 

We can remark here that from the Q it follows that {ui,U2)^^ ~ 
Uk and as we know one of its coefficient, we can compute it from the Uk- 
Analyzing the view of matrices Mj, 5*2"^ < i < S2 — I (namely, the presence 
of zero's on the "leading" or "trailing" diagonals) we see that we can fix the 
columns in such a way that {ui,U2y = for that i, so the structure of the 
sequence of is analogue to the one of usual subresultants. 

From the (j21) it follows the law of hk transformation: 

i^k+2 — iT'k+igk-i'^k+idk — i^k+igk-i-j— xTTfi^fc 

hkXi 9k-i gk'Si 



From the considerations above we can derive the algorithms for computing 
the gcd and resultants. Bellow we present the algorithm for gcd computa- 
tion, (we present in the style a la Algorithm C from ^Ij): 

Algorithm C 

Input: f,g are polynomials 
Output: the gcd of / and g 

C'l. [Reduce to full and primitive.] (u,v): = (f,g), d := gcd(cont(u), cont(v)), 
e := min(trailDeg^('u), trailDeg^(f )), replace {u,v) by 

(primpart(M)/a;*'''*''°°s^*^"\pr impart (f)/x*'"'^'^^'^s^('')). Ifdeg3,(M) < deg2,(v) then 
replace {u,v) by iv,u). Set h := 1, g := 1, g = 1, G := 1, G := 1. 
C'2. [General pseudo-remainder.] Apply genPRem(u, v) and assign r, 6, A, 
92, g2, w. If r = 0, then return (ix^f/cont(v). 

C'3. [Adjust remainder.] u := v; v := {rG) / {Ggh^) g := g2] g := g2] 
h := Gg^/{Gh^-^)- G := g^; G := g\ go to C'2 
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In the algorithm for computing the resultant of two full polynomials the 
algorithm is almost the same, but one need to return the value of h. 

For the non-full polynomials the following formula for resultant can be 
used: YeSx{xu,v) = tCx{v)TeSx{u,v) (up to the sign, of course). 

3 Implementation 

The algorithms for gcd and resultant computing above was implemented 
with the Axiomxl computer algebra system, which allows to get an efficient 
executing code. As a coefficient ring it was used the ring of polynomials 
Z[?/]. In Axiomxl there are two different structures for dense and sparse 
polynomials. As a relativeSize it was used the degree for dense polynomials 
and number of non-zero terms for sparse polynomials. The results of testing 
is the following: in the case of dense polynomials the algorithm is not slower 
than the usual subresultant algorithm; on some examples it is times faster 
than the usual subresltant algorithm. 

4 One property of generalized subresultants. 

Here we change the notation and will notate the generalized subresultants as 
SI to underline the analogues with usual subresultants. SI means that we get 
the generalized subresultant from the matrix for the usual subresultant Sk- 
The well known property of usual subresultants is that theire formal leading 
coefficients {principal resultants) flCx{Sk) allows one to check the degree of 
gcd [2|. The generalized subresultants have the same property, namely, the 
following lemma can be proved: 

Lemma. Let 5*^ be the sequence of generalized subresultants of two full 
polynomials A and B. Then deg^.(gcd(A, 5)) = d iS (flcx.(^^) or ftc^(^^)) = 
■ • ■ = (flCxiSU) or itCx{Sd-i)) = and flc,(S^) ^ (then also ftc,(^^) ^ 
and back); here flc^; and ftc^; are formal leading and trailing coefficient, they 
are some determinants. 

The content and proof of the lemma is almost analogous to the corollary 
7.7.9 from [2]. We just make here some remarks. Everywhere in the previous 
to the corollary 7.7.9 lemma's in [2| the PSCj appears it can be substituted 
by formal leading or trailing coefficient of the generalized subresultants. The 
big role in the proof plays the equation A{x)Tj{x) + B{x)Uj{x) = Cj{x), 
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where there is some conditions on the degrees of Tj{x), Uj{x) and Cj{x) and 
which holds when formal leading coefficient of Sj is vanishes. In the our case 
this equation will be of the form A{x)Tj{x) + B[x)Uj{x) = x*Cj{x), where 
X* means some power of x. 
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